Steps in the bacterial flagellar motor 
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The bacterial flagellar motor is a highly efficient rotary machine used by many bacteria to propel 
themselves. It has recently been shown that at low speeds its rotation proceeds in steps [Sowa 
et al. (2005) Nature 437, 916-919]. Here we propose a simple physical model that accounts for 
this stepping behavior as a random walk in a tilted corrugated potential that combines torque and 
contact forces. We argue that the absolute angular position of the rotor is crucial for understanding 
step properties, and show this hypothesis to be consistent with the available data, in particular the 
observation that backward steps are smaller on average than forward steps. Our model also predicts 
a sublinear torque-speed relationship at low torque, and a peak in rotor diffusion as a function of 
torque. 



Bacteria swim by virtue of tiny rotary motors that 
drive rotation of helical fiagella. These motors are pow- 
ered by a transmembrane proton (or Na+) flux which is 
converted into torque. However, little is known about 
the detailed mechanisms of energy conversion, or torque 
generation. Recently, a new result has provided direct 
insight into motor operation [1]: at low speeds, the bac- 
terial flagellar motor proceeds by steps. This stepping 
is stochastic in nature, as manifested by the occurrence 
of occasional backward steps even for motors locked in 
one rotation direction. What is the origin of motor steps 
and how can these steps be reconciled with the near per- 
fect efficiency of the motor observed at low speeds [2]? 
We argue that steps, including backward steps, are an 
inevitable consequence of the physical structure of the 
motor — a stator driving a "bumpy" rotor through a vis- 
cous medium. 

In response to chemotactic signals, flagellar motors 
switch from counterclockwise to clockwise rotation caus- 
ing cells to tumble or change directions. In E. coli, 
the basic mechanism of torque generation appears to 
be the same for both directions of motor rotation [3]. 
Torque is generated by the passage of H+ ions (or in 
some organisms Na+ ions) through the cytoplasmic mem- 
brane. As shown schematically in Fig. lA, torque is ap- 
plied to the rotor, including the fiagellum, by the sta- 
tor, which is comprised of independent torque-generating 
units (MotA/B complexes) anchored to the peptidogly- 
can cell wall. The exact number of torque- generating 
units can vary from motor to motor, with the maxi- 
mum estimated to be at least 11 [4]. The rotor in- 



cludes 26 circularly arrayed FliG proteins that contact 
the MotA/B complexes. The torque-speed relation of 
the motor has been measured under a range of conditions 
[5, 6, 7, 8]. The maximum torque in the high load, low 
speed regime tracks the electrochemical potential differ- 
ence or proton motive force (PMF) across the membrane, 
and the motor operates with nearly perfect efficiency [2] . 
Whereas torque and efficiency fall off at high speeds, pro- 
ton flux and motor rotation are always strongly coupled 
with ^ 120 protons passing through the membrane per 
MotA/B unit per rotation [2]. 

Recent experiments, where rotation was measured 
by attaching a polystyrene bead to a flagellar stump 
driven by a counterclockwise-locked Na+-powered chi- 
maeric motor at low speeds (low Na+ -electrochemical- 
potential difference and low stator number), revealed 
that the motor proceeds by steps [1]. The steps have 
average size ^ 13.8°, which corresponds to 26 steps per 
rotation, exactly the number of copies of FliG around 
the rotor. Occasional backward steps are observed and, 
interestingly, these are smaller on average than forward 
steps (10.9° versus 13.8°). These observations, as well as 
the stepping mechanism itself, have so far remained un- 
explained. It has been suggested that stepping is caused 
by the stochastic passage of ions. However, as pointed 
out in [1], the energy provided by passage of a single ion 
can only move the rotor attached to a 1 /im polystyrene 
bead by 5°, much less than the typical observed step size. 

Here we propose a simple physical model to explain 
stepping: the stator applies nearly constant torque to the 
rotor, but, at the same time, contact forces on the rotor 
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produce a potential and therefore an additional torque 
with approximately the 26- fold periodicity of FliG. Our 
model naturally accounts for the existence of backward 
steps, as well as the discrepancy between forward and 
backward step sizes, and also predicts that step statistics 
depend on the absolute position of the rotor around the 
circle. Our predictions are found to be consistent with 
the available data, including angular diffusion of the mo- 
tor [9], and suggest how steps could be used to study 
the physical structure of the motor. A novel testable 
prediction is that the torque-speed relation will become 
sublinear at very low torques. 



Model 

Our model for stepping relies on two main assump- 
tions: constant or nearly constant torque between stator 
and rotor and an approximately 26-fold periodic contact 
potential. 

(1) All torque- generating units apply torque simultane- 
ously and additively. Following the model of Meister et 
al [10], we assume that each MotA/B complex acts as 
a set of protein springs that reversibly store the energy 
available from H+ (or Na+) translocations (see Fig. IB). 
The protein springs are attached to fixed sites of the 
rotor circumference. When an H+ passes through the 
membrane, it causes a spring to detach from its attach- 
ment site, stretch, and reattach to the next site. At 
stall, all springs are maximally stretched, such that the 
PMF matches the energy necessary to stretch a spring 
to its next site. At low speeds, the rotor moves and 
springs relax, but these are quickly restretched by H+ 
passage, so that the system remains in quasi- equilibrium 
with the torque set by the PMF. Spring stretching may 
vary slightly among units, but since there are several mo- 
tor units, the instantaneous torque self- averages and is 
nearly constant in time. Under this scenario, steps can- 
not be explained at the level of a single pair of MotA/B 
and FliG subunits, but must arise at the global level of 
the rotor-stator interaction. 

(2) There are contact forces between the stator and the 
rotor. These forces may be caused by contact between 
the MotA/B stator units and FliG proteins, but also pos- 
sibly by contact with FliF (M-S ring), FlgG, FlgH or Flgl 
proteins (distal rod, L, and P rings, respectively) each of 
which forms a circle of 26 copies. There may be other 
periodicities to the contact forces as well, arising from 
the filament and the hook, which are 11-fold periodic, 
from FlgK and FlgL (hook- filament junction, 11 copies 
each), FlgB, FlgC and FlgF (proximal rod between L 
and P rings, 6 copies each) and FliE (Rod-MS-ring junc- 
tion, 9 copies). We assume that a 26- fold periodicity is 
dominant, in agreement with experimental observations. 
We therefore collect all contact forces in a potential V{0) 



which we suppose to be approximately 26-fold periodic 
(Figure IC). 

Since the motor operates at the molecular scale, its 
rotation is intrinsically stochastic as it is subject to ran- 
dom thermal fluctuations. Under the combined influence 
of the applied torque, the contact potential, and ther- 
mal fluctuations the rotor performs a continuous random 
walk in a tilted, approximately periodic potential, which 
we model by the following Langevin equation: 



dO 
dt 



IdU 



(1) 



where v is the drag coefficient, r the total torque exerted 
by the stator, and where the potential U {0) includes both 
the torque and the contact potential: 



U{0)= V{0) - tO. 

(contact potential) (torque) 



(2) 



The term ^(t) represents Gaussian white noise and ac- 
counts for thermal fluctuations: (C(^)C(^O) = 2DS{t — f), 
where D is the rotor diffusion coefficient, related to the 
temperature and the drag coefficient via Einstein's re- 
lation: D = kT/v. In experiments, a load (usually a 
polystyrene bead) is attached to a flagellar stump and 
this load is largely responsible for the drag. For simplic- 
ity, we assume that linkage between motor and load is 
instantaneous, as the relaxation is rapid compared to the 
typical stepping time (see Discussion). 



Results 



Steps are barrier-crossing events. Numerical simu- 
lation of the model (Eq. 1) shows that rotation proceeds 
in steps (Fig. IC). These steps correspond to jumps be- 
tween adjacent wells of the tilted potential. Jumps/steps 
are possible thanks to thermal fluctuations which drive 
the system out of energy minima; without these fluctua- 
tions, the system would remain stuck in one well forever. 
Steps therefore correspond to crossings of the energy bar- 
riers separating wells. According to the Arrhenius law, 
the average time to cross a barrier increases exponen- 
tially with the barrier height. Because of the tilt induced 
by the torque, steps in the forward direction correspond 
to lower energy barriers than steps in the backward di- 
rection (c/. Fig. IC). Forward steps are therefore more 
likely to occur than backward steps, so that on average 
the motor moves forward. 

Backward steps are smaller than forward steps. 

To further investigate stepping in our model, we wrote a 
step detector algorithm similar to that described in [1], 
and applied it to a simulation of ?^ 2 • 10^ rotations (see 
Methods). Fig. 2 A shows the histogram of step sizes. 
As in the experiment, we flnd that backward steps are 
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FIG. 1: Model for stepping of the flagellar motor. A. Side view of the flagellar motor. B. Top view of the motor 
highlighting the model's essential ingredients. The passage of across the inner membrane causes the stretching of protein 
"springs" which link the fixed stator complexes (MotA/B) to the rotor (FliG, etc.). The stretched springs apply a torque to 
the rotor. Contact forces between the stators and the rotor also produce a potential of interaction, which is approximately 
26- fold periodic due to the 26 FliG subunits. C. Left: Rotation of the rotor as a whole corresponds to a viscously damped 
random walk in a tilted corrugated potential IJ(&) arising from the combined torque and contact potential. Right: Example of 
a trace generated by the model (blue) and the inferred steps (red) between local potential wells (shown with purple shading). 



smaller on average than forward steps: the mean forward 
step is ^ 13.8° 360°/26), against 12.0° for backward 
steps. (The precise values vary with the particular choice 
of potential and torque, but the mean step size is always 
larger for forward than backward steps.) 

Recognizing steps as barrier crossing events allows us 
to readily explain the difference between average forward 
and backward steps sizes using a simple intuitive argu- 
ment, as illustrated in the inset to Fig. 2 A. Backward 
steps occur infrequently because the energy barriers for 
these steps are higher than for forward steps. By con- 
trast, all forward steps must occur as the motor moves 
forward, regardless of the heights of energy barriers. In 
addition, as shown in the inset to Fig. 2A, barrier heights 
and step sizes tend to be positively correlated. Roughly 
speaking, higher barriers extend over longer ranges. (We 
confirmed the generality of this correlation by simulating 
motion for many approximately 26-fold periodic poten- 
tials V(S)^ with the results shown in Fig. SI.) Therefore, 
backward steps occur mostly over the lower barriers, and 
lower barriers correspond to smaller step sizes. This im- 
plies that backward steps are on average smaller than 
forward steps. Note that this argument relies on the fact 
that the barriers are not all identical. 

To test the scenario proposed above, we asked whether 
the size of forward steps immediately preceding or fol- 
lowing backward steps differ from the average of 13.8°. 
According to our picture, the barrier crossed by these for- 
ward steps should be the same as the one crossed by the 
backward step immediately preceding or following, im- 
plying a small barrier and therefore a small forward step 
size. Fig. 2B-E shows that indeed forward steps preced- 
ing of following a backward step are smaller on average, 
12.0° in the experiment, 12.4° in our simulation, than 



the mean forward step of 13.8°. (Note that even for for- 
ward steps over the same barrier, backward steps are still 
slightly smaller on average in both experiment and sim- 
ulation. This suggests that the step detection algorithm 
has a small systematic bias — see supporting information 
(SI) Text.) 

Importance of absolute position of the rotor. The 

difference between forward and backward step sizes relies 
on the 26 barriers around the circle not all being iden- 
tical. These heterogeneities may exist because the po- 
tential V(0) is only approximately 26-fold periodic and 
contains other periodicities as well arising from the fila- 
ment, hook, or other parts of the rotor. In any event, an 
essential prediction of our model is that step sizes and 
backward step frequencies will depend strongly on ab- 
solute position {i.e. modulo 360°), reflecting the fixed 
contact potential V(S). We now examine how step fre- 
quencies and sizes depend on the properties of particular 
barriers, specified by the position of the rotor around 
the circle, and how this can tell us something about the 
detailed nature of contact forces. 

According to our model, backward steps should be 
much more likely to occur at low barriers. In contrast, 
where forward steps occur should be much less sensitive 
to barrier heights. This follows simply because for each 
complete rotation the number of forward steps over any 
barrier is one plus the number of backward steps over 
that same barrier. Therefore, as long as backward steps 
are rare, the average number of forward steps over each 
barrier will be close to one per rotation and therefore the 
frequencies of forward steps will be similar for all barri- 
ers. This is illustrated in Fig. 3A, which presents average 
step frequencies for each of the 26 barriers of a particu- 
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FIG. 2: Backward steps are smaller than forward 
steps. A. Probability distribution of step sizes using a ap- 
proximately 26-fold periodic potential (see main text), show- 
ing that backward steps are on average smaller than forward 
steps, in agreement with experiment. Inset: Backward steps 
rely on low barriers, which occur preferentially where angular 
steps sizes (e.^. A^2) are small. B-E. Forward steps imme- 
diately following or preceding backward steps are found to 
be smaller on average (12.0° in the experiment, 12.4° in the 
model) than the mean of all forward steps (13.8° in model and 
experiment, black dots). Crosses denote mean and standard 
deviation of backward and subsequent or previous forward 
steps, while black dots and horizontal lines give the mean 
(13.8°) and standard deviations of all forward steps. Note 
change of scale between simulation and experiment. The ex- 
perimental data are from [1]. 



lar potential (chosen to be the same as in Fig. 2), as 
calculated from first-passage theory (see Methods). As 
expected, there is considerable variation among barriers 
in backward-step frequencies, but much less variation in 
forward- step frequencies. 

To relate average step sizes to absolute position, we 
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FIG. 3: Absolute position of the rotor matters. A. A 

typical approximately 26-fold periodic tilted potential IJ (red 
line) used in our simulations. For each of the 26 barriers in 
[/, we show the frequencies of forward and backward steps 
across that barrier obtained analytically from first-passage 
theory (see Methods). Backward steps occur much more fre- 
quently at low barriers. B. Average backward and forward 
step sizes for each of the 26 barriers around the circle. Each 
color corresponds to a simulation with a different choice of 
the potential IJ . C. The ratio of backward over forward step 
counts for a given barrier decreases with the average forward 
step size (colors as in B). D and E. Same as B and C, but 
with experimental data [1]. Each color corresponds to a dif- 
ferent cell. For each cell and each position around the circle, 
we show a data point only if there were at least 10 backward 
steps. 



examined the sizes of backward and forward steps for 
each of the 26 barriers. To properly assign steps to bar- 
riers, we sorted steps into 26 equal bins according to the 
angular position of the rotor when the step occurred, 
and calculated the average backward and forward step 
sizes in each bin. We applied this procedure to simula- 
tions of ~ 10^ rotations generated with three different 
potentials (see Methods), and to four experimental 
traces, corresponding to four distinct cells, totaling 700 
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rotations [1]. According to our model, forward and back- 
ward steps across the same barrier should have the same 
average size. In the simulations (Fig. 3B), mean forward 
and backward steps across the same barrier were found 
to be linearly correlated, though with a systematic offset 
toward smaller backward steps. (As discussed above, this 
offset is likely the result of a bias in the step detection 
algorithm.) We found that mean forward and backward 
step sizes across the same barrier are also positively corre- 
lated in the experiment (Fig. 3D), in agreement with our 
prediction (with the same bias towards smaller backward 
steps). Overall, these results suggest that the absolute 
position of the rotor accounts for much of the variability 
observed in step size, and supports our model of a fixed, 
nearly periodic contact potential. 

We next show that, both in the simulation and in the 
experiment, the barriers with a high frequency of back- 
ward steps are the same barriers where step sizes are 
short. To this end we plot the frequency ratio of back- 
ward steps to forward steps across each barrier versus the 
average forward step size (which we showed in Fig. 3B,D 
correlates with the backward step size), for both simula- 
tions (Fig. 3C) and experiment (Fig. 3E). In both cases 
backward-step frequencies fall off sharply with average 
forward step size. 

Since barriers where steps are smaller have higher 
backward-step frequencies than other barriers, they con- 
tribute more to the average backward step size. There- 
fore, the mean backward step is smaller than 360°/26 ~ 
13.8°, which would be the mean backward step size if all 
barriers contributed equally. In contrast, forward- step 
frequencies vary little from barrier to barrier, so that 
all barriers contribute more or less equally to the aver- 
age forward step size, which is therefore approximately 
13.8°. This explains why backward steps are smaller than 
forward steps on average. 

Sublinear torque-speed relation. The recognition 
that steps are barrier-crossing events has a direct im- 
plication for how rotation speed depends on torque. In 
the absence of contact forces {i.e. V{0) = 0), the average 
rotation speed / depends linearly on torque: / = r jl-KV. 
However, when the contact forces are comparable to the 
torque, rotation is hindered by barriers, and the system 
spends much of the time in local energy minima. Ro- 
tation is then not only limited by drag, but also by the 
rate of barrier crossing, leading to lower rotation speeds: 
/ < r/27ri/. We computed analytically (see Methods) the 
torque-speed relation for loads with various drag coeffi- 
cients for a perfectly 26-fold periodic sinusoidal potential 
with amplitude 1.5A:T, as shown in Fig. 4A (using approx- 
imately 26-fold periodic potentials yielded qualitatively 
identical results). At high torques, the linear relation is 
recovered asymptotically, which follows because increas- 
ing torque decreases the barriers to forward rotation, and 
eventually eliminates them completely, as shown in the 



insets to Fig. 4A. 

Rotor diffusion. Fig. 4B shows the effect of the con- 
tact potential on the effective long-time diffusion coef- 
ficient Deff- At low torques, diffusion is slowed down 
by barriers, while at high torques one recovers the natu- 
ral diffusion coefficient D. Interestingly, at intermediate 
torques rotor diffusion is actually enhanced by the con- 
tact potential. In this regime, the contact potential is a 
small but variable correction to the torque. This vari- 
ability contributes to the variance of the rotation speed, 
thus effectively enhancing rotor diffusion (see SI Text). 
At large torques, Deff approaches D asymptotically: 



1 + 3 



{{deVf)e 



(3) 



To check the consistency of the predicted diffusive be- 
havior against previous results, we compared our model's 
prediction with experimental measurements of the vari- 
ance in the rotation time [9, 11]. In [9], the rotation time 
of a tethered cell was measured. Simple diffusion predicts 
that the variance in rotation time per cycle is: 
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For a single fully-powered torque- generating unit, torque 
is estimated to be r = 250 pN.nm.rad"-^ [7]. The 
measured speed in [9] was noisy and depended on the 
particular cell but was about / 3 — 5 Hz for three 
torque-generating units, leading to an estimate for the 
drag coefficient v = r/(27r/) 25 — 40 pN.nm.s.rad"^, 
and for the diffusion constant D = kT/v « 0.1 — 0.15 
rad^/s. Thus the variance in cycle time for three torque- 
generating units predicted by simple diffusion is (^T^) ^ 
0.7 — 2.0 • 10~^ s^, which is consistent with the reported 
value of ((5T2) ?^ 2.1 • 10"^ s^ [9]. We conclude that it 
may not be necessary to consider other sources of fluctua- 
tions {e.g. proton translocations) to explain the observed 
variance in cycle time. Note that in these experiments 
the torque was high, and therefore contact forces are not 
expected to have had a significant effect on diffusion. 

Our model also predicts a negative feedback reaction 
from the MotA/B protein springs that in principle could 
reduce diffusion. Namely, every time the rotor moves 
forward the springs relax causing a transient decrease of 
torque, and the opposite every time the motor moves 
backward. To estimate the magnitude of this effect, we 
model torque dynamics by linking spring elongation to 
rotor position, and by assuming that the springs "re- 
stretch" to their equilibrium position prescribed by the 
PMF with a characteristic relaxation time (see Meth- 
ods). Within this model, the effective diffusion coeffi- 
cient is found to be D{1 — 2/i/z/), where ji is the slope 
of the torque-speed relation of the motor near stall. The 
value of /i ranges from .03 to .07 pN.nm.s.rad"^ depend- 
ing on the temperature [6] and is therefore much smaller 
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FIG. 4: The model predicts a sublinear torque-speed relation and a peak in rotor diffusion. A. Rotation speed of 
the rotor as a function of torque for loads with different drag taken from [7]: solid curves, from top to bottom, 27viy = 2.7, 4.3, 
5.1, 8.5, and 28 pN.nm.s.rad"^. The potential was chosen to be perfectly 26-fold periodic: V(0) — Asin(26^), with A — 1.5 /cT. 
As torque increases, the rotation speed asymptotes to the behavior expected in the absence of barriers, / = T/27rz/, represented 
here for 27rz/ = 2.7 pN.nm.s.rad"^ (dashed line). At small torques, the rotation speed is limited by the rate of barrier crossing 
(left inset), while at high torques the tilt makes barriers easy to cross (right inset), and rotation is only limited by drag. B. 
Effective diffusion coefficient as a function of torque for a load with drag coefficient 27rz/ = 1 pN.nm.s.rad"^: solid curves, 
V{0) = Asin(26^) with A = 0.5 kT, 1 kT, and 1.5 kT. The real diffusion coefficient D is represented by a dashed line. 



than relevant values of the drag coefficient u ^ 0.5 — 50 
pN.nm.s.rad"^. We conclude that in the conditions of 
the discussed experiments, the effect of negative feedback 
from the springs on diffusion is negligible. 

Distribution of waiting times. According to our 
model, the distribution of waiting times between steps is 
expected to be roughly exponential. The distribution of 
waiting times between steps (forward or backward) for an 
exactly 26-fold periodic potential is indeed exponential 
(Fig. S2A). When the potential is heterogeneous, the av- 
erage waiting time depends on the barrier. Even though 
the distribution of waiting times across each barrier is 
exponential, the overall waiting-time distribution is not, 
appearing rather as a "stretched" exponential (Fig. S2B). 
The experimental distribution also resembles a stretched 
exponential (Fig. S2C). 

Discussion 

Our model explains stepping of the bacterial flagellar 
motor by interpreting its rotation as a viscously damped 
random walk driven by a constant torque and by a het- 
erogeneous contact potential caused by the physical irreg- 
ularities of the rotor. In this picture, steps are recognized 
as barrier-crossing events between adjacent minima of a 
tilted and corrugated energy potential. Corrugations are 
caused by contact between the stators and the protein ar- 



rays (FliG, among others) making up the rotor structure. 
Recently a more accurate picture of this structure has 
emerged, thanks notably to electron microscopy studies 
[12]. 

Our model predicts a 27r periodicity of the potential, 
so that the absolute angular position of the rotor with 
respect to the stator is an underlying determinant for 
step statistics, and this prediction is found to be consis- 
tent with the available experimental data. In particular, 
our model offers an explanation for the experimental ob- 
servation that backward steps are smaller than forward 
steps on average. 

Another prediction of the model is that rotor speed 
grows sublinear ly as a function of torque. At low torques 
rotation is slow because of trapping in local minima, 
whereas at high torques the barriers between minima are 
lowered and eventually eliminated. Additionally, we pre- 
dict that at low torques rotor diffusion is hindered by 
barriers, while at high torques the variability of the po- 
tential actually enhances diffusion. Although in principle 
other sources of fluctuations, such as ion translocation, 
could impact rotor diffusion, we showed that in the rel- 
evant regimes simple diffusion can account for nearly all 
of the observed variance in cycle time [9]. In order to 
verify these predictions experimentally, one would need 
to simultaneously measure the rotor speed and the pro- 
ton (or Na+) motive force, believed to be proportional 
to torque at low speeds, in the regime where torque and 
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contact forces are comparable. (Note that in the step- 
ping data [1] we have analyzed, torque could vary during 
the course of the experiment as the result of changes in 
the number of motor units.) 

Our model is consistent with other experimental re- 
sults on the bacterial flagellar motor. Because the model 
relies on the assumption that the energy from ion translo- 
cation is reversibly stored in protein springs [10], it im- 
plies a near-perfect efficiency of the motor at low torques 
[2]. The same mechanism can account for both clock- 
wise and counterclockwise motor rotation [3] — these two 
cases simply corresponding to the springs being stretched 
in opposite directions. If the contact potential stays the 
same when the motor changes direction, our model pre- 
dicts that backward steps will occur preferentially at the 
same absolute angles irrespectively of the direction of ro- 
tation. The observation that the duty ratio is very close 
to one even with a single torque-generating unit [7] can be 
encompassed in our model by assuming that each torque- 
generating unit comprises at least two springs. 

In our analysis we have neglected one effect that is 
not crucial for our analysis, but which may prove impor- 
tant for inferring the detailed nature of the contact po- 
tential. Specifically, we have assumed that equilibration 
of the elastic linkage between the motor and the load 
is rapid compared to the waiting time between steps. 
For a torsion constant ke ^ 400 pN.nm.rad"^ [13, 14] 
between the rotor and the load, and a drag coefficient 
27ru ~ 1.0 — 10 pN.nm.s.rad"^, the relaxation time is 
te = ly/ko = 0.4 — 4 ms. In contrast, the typical waiting 
time between steps ranges from 10 to 100 ms, depending 
on experimental conditions. If the elastic linkage was too 
soft, the polystyrene bead would respond to the motion 
of the rotor with a delay te , and steps would be smoothed 
out. This does not seem to occur in the experiment. 

Another effect, which we have considered (see Rotor 
diffusion) but did not include in our simulations, is the 
relaxation of MotA/B protein springs as a rotor step 
occurs. For example, when the rotor moves forward, 
the torque decreases because the protein springs relax. 
Usually these springs are restretched so quickly by ion 
translocation that the transient decrease of torque can be 
neglected. However, during a barrier crossing event the 
rotor motion might be so fast that protons are not able 
keep up. This would result in a temporary drop in torque 
and make barrier crossing more difficult. A similar argu- 
ment applies to backward steps. We have already shown 
that at the "mean-field" level, where rotation speed and 
ion flux are time- averaged, this negative feedback has 
only a small effect. However, the instantaneous rotation 
speed during a step can be much larger than its average. 
How fast can a proton translocate through the motor? 
The maximum flux of protons through a single motor 
unit can be estimated by considering the maximum rota- 
tion speed before the torque starts dropping (the "knee" 



of the torque-speed relationship [5]). For a single motor 
unit in natural conditions, this speed is about 150 Hz for 
a torque of 250 pN.nm.rad"^ [7]. The power generated by 
the motor is then 150 x 27rx 250 ^ 240,000 pN.nm.s-^ 
Each proton provides at most 150 mV xe = 24 pN.nm, 
so that the number of protons per second is at least 
240,000/24 = 10,000. The timescale of proton passage 
is therefore less than 0.1 ms. A single rotor step corre- 
sponds to the passage of 3 protons [1], so restretching 
the protein springs should take less than 0.3 ms, which 
is below the current experimental time resolution. For 
comparison, an instanton calculation [15] reveals that 
the typical time for crossing a barrier is bounded from 
below by (27r/26) x z//max|S^C/| ^ 0.2 ms (for 27rz/ = 1 
pN.nm.s.rad-\ and U{0)/kT = 1.5cos(266>) - 100). 

Other molecular motors have shown stepping behav- 
ior, including the actin- myosin motor [16], the dynein- 
microtubule motor [17], and kinesin [18]. In these ATP- 
powered motors, which are less powerful than the bacte- 
rial flagellar motor by orders of magnitude, stepping is a 
built-in and essential part of motor operation. By con- 
trast, we have argued that in the bacterial flagellar motor 
the observed stepping arises solely from steric hindrance. 

Our work leaves open a number of questions. It would 
be interesting to infer the precise form of the contact po- 
tential V{0) from rotation data and see how and whether 
it varies in time and among motors. To this end a more 
sophisticated approach to learning the potential may be 
required, e.g. employing maximum likelihood techniques. 
Lastly, one still needs to understand the mechanism of 
torque generation, including the role played by the dis- 
creteness of ion translocation, the chemical nature of pro- 
tein springs and their attachment sites, as well as the 
energy conversion process. 



Methods 

All the experimental data presented in this paper is 
from [1] and were used with the kind permission of 
Richard Berry. 

The simulation data were obtained by numerical in- 
tegration of Eq. 1 by Euler's method. The same step 
finding algorithm as the one described in [1] was used on 
both simulation data and experimental data to extract 
steps. 

The step frequencies presented in Fig. 3 A were ob- 
tained analytically using ffist-passage theory [19, 20]: La- 
bel the wells i = 1, . . . , 26, and denote the local minima 
of U{0) by Oi. Consider three consecutive wells centered 
at Oi-i, Oi and respectively. Starting at 0^, call 
the probability of first jumping forward and 7r~ the prob- 
ability of first jumping backward. These probabilities are 
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given by: 



mean- field differential equations: 



,u(e)/kT 



(5) 



Given the transition probabilities tt^ and 7r^~, we write 
a master equation for the probabihty pi (m) of the rotor 
being in weU i after k steps. 



Pi{m) =pi_i{m- l)7r^ti Pi+i{m - l)7r. 



(6) 



At large m a stationary state is reached, and pi := 
Pi{m oo) satisfies the conservation equation: 



Pi =Pi-lT^t-l ^Pi+l'^i+V 



(7) 



The step frequencies presented in Fig. 3 are then given 
by piTi^ and Piii^ . (Note that the step frequencies sum 

to one, Ei=iPi< + Y.i=iPi^i = !•) 

The torque-speed relation shown in Fig. 4 A was also 
estimated using first-passage theory. For simplicity we 
assumed a perfectly 26-fold periodic potential, but the 
results are qualitatively the same when the periodicity is 
only approximate. The average time to move from one 
minimum of the potential to the next is given by: 

(4tep) - 1 _^ g-rA//cT ' 

where A = 27r/26 and r is the torque. This yields a 
rotation speed: 



/ 



A(7r+ 



27r(istep) 



— ; T tanh — — ; 

27r(4tep) \2kT 



(9) 



where 7r+ and tt are obtained from Eq. 5. The asymp- 
totic expansion of Eq. 9 for r \dV/dO\ , kT yields 



T 

2m' 



(10) 



The effective diffusion coefficient Dgfr shown in Fig. 4B 
is estimated using similar techniques (c/. [19, 20] and SI 
Text). 

To estimate the magnitude of the negative feedback 
of the MotA/B protein springs on diffusion, we model 
the dynamics of rotor angle and torque by the following 



de 

dt 
dr 

'dt 



' ^PMF 



^relax 



k 



de 

di' 



(11) 
(12) 



The first equation is the same as Eq. 1 but with a vari- 
able torque and without contact forces. In the sec- 
ond equation, torque is assumed to follow the stretch- 
ing/unstretching of the springs with rotation, and there- 
fore the rate of change of torque is linearly related to the 
rate of change of rotor angle through an effective spring 
constant k (second term of r.h.s. of Eq. 12). At the 
same time, due to the restretching of springs upon pro- 
ton passage, torque relaxes toward its equilibrium value 
TpMF (first term of r.h.s. of Eq. 12). Solving the sec- 
ond equation for steady-state rotation yields the torque- 
speed relationship for the motor: r = Tpmf — 27r/i/, 
with /i := /ctreiax- Solviug both equations for the ef- 
fective diffusion coefficient in the limit /i/i/ <C 1, we find 
{Se^)/t = 2D [1 - 2fi/u + Oiifi/u)^)] . 

For all numerical simulations and analytic calculations 
the potential was chosen to be of the form: 

V{0) = ylcos(26(9) -h5cos(10(9) -hCcos(ll(9). (13) 

Experimentally, there is evidence for components of the 
contact potential with ~ 10 — 11 fold periodicity (see 
Ref. [1], Fig. 3b). Except where stated otherwise, we 
used A = 1.5A:T, T = 290Kand 27ri/ = 1 pN.nm.s.rad-^ 
In Fig. 2 and 3A, we used B = C = 0.6 kT. In Fig. 3B 
and 3C we chose three sets of values for C, and r, 
in units of kT: {B = 0.6, C = 0.6, r = 10), (5 = 
1, C = 0.5, T = 15) and {B = 0.3, C = 0.8, r = 12). In 
Fig. 4, the torque-speed relation and rotor diffusion were 
calculated with B = C = 0. 
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Supporting Information 



Step finding algorithm. Steps were found using the 
algorithm described in [1]. First, an entire episode (angu- 
lar position vs. time) was fitted by a single step function, 
and thus divided into two intervals. This procedure was 
repeated iteratively iVgteps times — at each iteration the 
interval for which the data had the largest range of an- 
gles was replaced by a best-fit step function. Then, a 
quality factor Q = {x2 — xi)^/(var xi/ni + var X2/n2) 
was calculated for each assigned step, where xi (^2) is 
the angular position left (right) of the step, and ni (712) 
is the number of data points on the left (right) of the 
step. Finally, steps with a low quality factor were re- 
moved and their adjacent intervals merged until all steps 
have a quality factor greater than Qmin- 

This step finding algorithm introduces a small bias that 
tends to underestimate the sizes of backward steps. For 
each interval, the algorithm finds the mean value of the 
angular position in this interval. But steps do not occur 
instantaneously, and the data points leading from one in- 
terval to the other are themselves included in the interval 
means. As a result, these means are biased toward where 
the rotor is coming from and where it is going. When the 
rotor is stepping forward, these two biases tend to can- 
cel each other. However, when the rotor steps backward, 
these steps are usually both preceded and followed by 
forward steps. Thus, the interval before the backward 
step is biased toward lower angular positions, and the 
interval after the backward step is biased toward higher 
positions. As a result, the biases reinforce each other 
such that backward steps are estimated to be smaller 
than forward steps. We have checked that even on a 
perfectly 26- fold periodic potential V{0) = 74cos(26^), 
with A = l.bkT and a torque r = 10/cT, the most prob- 
able forward step size according to the algorithm was 
13.8° ^ 360°/ 26, which is correct, while the most proba- 
ble backward step size was 13.0°, which is 0.8° to small. 

Effective diffusion. The effective diffusion coefficient is 
estimated using the techniques of first-passage theory [2] . 
For an exactly 26-fold periodic potential, this effective 
diffusion coefficient is given by [3]: 



.2 r 



1 



(tH 

(^step) (^1 



+ 



step 



step 



> - 2(tstep>') 



(1) 

where A = 27r/26 is the step size. 7r+ and 7r~ are the 



probabilities of first jumping forward or backward, re- 
spectively, and are given by: 

^+ - ^- = tanh j , ^+ + = 1. (2) 

tstep is the time before a step occurs (forward or back- 
ward). Its first two moments are: 



(^step) 



1 d0 e-^W/fe^ de' e^(^')/fe^ 
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)-u(e")]/kT 



(/.'^ d0' e^^C'')/'^^) d0' de" emo')-u(0")]/kT 



D J_^d0' eU(.e')/kT 



Expanding Eq. 1 for large torques gives: 



1 + 3- 



(5) 



(6) 



Waiting-time distribution. We recorded the distribu- 
tion of waiting times obtained by the step-finding algo- 
rithm, both for the simulation and for the experiment. 
For the simulation, in the case of a perfectly 26-fold pe- 
riodic potential (Fig. S2A), there is only one type of bar- 
rier, and the waiting-time distribution is approximately 
exponential. When the potential is only approximately 
26- fold periodic (Fig. S2B), the waiting-time distribution 
is the sum of 26 exponentials, and resembles a stretched 
exponential. The waiting-time distribution for the ex- 
periment (one cell. Fig. S2C) is consistent with a sum 
of exponentials. The stretched appearance of the experi- 
mental distribution may be due to non- uniform barriers, 
as predicted by our model, but may also be due in part 
to the observed variability in average speed, presumably 
due to changes in torque, during the course of the exper- 
iment. 



[1] Yoshiyuki Sowa, Alexander D Rowe, Mark C Leake, Toshi- 
haru Yakushi, Michio Homma, Akihiko Ishijima, and 
Richard M Berry. Direct observation of steps in rotation of 
the bacterial flagellar motor. Nature, 437(7060):916-919, 
Oct 2005. 

[2] Crispin W. Gardiner. Handbook of Stochastic Methods: 

for Physics, Chemistry and the Natural Sciences (Springer 
Series in Synergetics). Springer, 3rd ed. edition, April 
2004. 

[3] B Lindner, M Kostur, and L Schimansky-Geier. Optimal 
diffusive transport in a tilted periodic potential. Fluct 
Noise Lett, 1(1):R25-R39, Jan 2001. 



3 



Q. 
Q) 
CO 

o 

03 
JD 

O 

g) 

03 
DQ 




12 13 14 15 

Step size (degrees) 



16 



17 



Fig. S 1: The height of barriers to backward steps is positively correlated with the backward step size. For 
100 randomly generated potentials V{0), we plotted the size of each of the 26 steps versus the height of the corresponding barrier 
to backward steps {i.e. the energy difference between the minimum of U{6) and its immediately preceding maximum). Each color 
corresponds to a particular potential. We used: V{0) = A26 cos(266') + Aq cos(66') + Ag cos(96>) + Aio cos(106') + An cos(116>) + 
AiQ cos(16^), with A26 — l.^kT, and Aq, . . . , Aiq drawn uniformly at random with fixed total power J2j — (1-0/cT)^. The torque 
is r = lO/cT ^ 41 pN.nm.rad"^ 



A. Simulation witli exactly 26-fold periodic potential 
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B. Simulation with approximately 26-fold periodic potential 
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Fig. S 2: The distribution of waiting times between steps (backward or forward). A. Simulation with a perfectly 
26-fold periodic contact potential V^(^). There is only one type of barrier, and the distribution of waiting times is roughly exponential. 
(The contact potential is V{0) = Asin(266') with A = 1.5kT, and torque r = 10A;T ^ 41 pN.nm.rad"^) B. Simulation with an 
approximately 26-fold periodic contact potential. There are 26 types of barriers, each of them having a different characteristic time. 
The overall waiting time distribution is a therefore the sum of 26 exponentials. (The contact potential and torque are the same as 
in Fig. 2 of the main text.) C. Experiment- the distribution is consistent with a sum of exponentials. 



